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SUMMARY 


Heat flow in a disc due to moving heat generated by electron 
beam was studied analytically. 

Computer programs based on the finite-element method were 
developed for the analysis of two- and three-dimensional mathematical 
models. 

The limited experimental data were compared with the analytical 
results and various factors which had influence on heat flow in the 
disc were studied. Among various factors radiation seemed to be 
important for the analysis of the disc which was stored in a chamber. 

The computer programs were modified to consider the effects of 
metal melting and solidification as well as of radiation. Temperature 
dependency of thermophysical properties was also considered in the 
program. 


SECTION I : INTRODUCTION 


The objective of the M551 experiment is to examine the molten 
flow and solidification of metal alloys during electron-beam melting 
in space. The specimens are circular discs of Stainless steel, 

2219 Aluminum and Tantalum, respectively. 

Of the various factors which had an effect on the metal melting 
and solidification in the experiment, gravity appears to have had 
the most significant effect, especially in space. However, gravity 
has no effect on heat conduction and radiation. 

The M.I.T, efforts placed emphasis on the heat flow analysis 
due to conduction and radiation. The information obtained by 
thermal analysis was intended tc assist other cortractors in their 
study on the metal melting and solidification. Mathematical models 
based on the finite element method were developed to analyze various 
heat-transfer phenomena involved in the M551 experiment. The com- 
puter program by use of the marhematical models was also developed 
to conduct heat flow analysis of a test specimen which had complex 
configurations and experimental conditions. 

Through the M551 experiment the M.I.T. efforts have been made 
in the following phases: 

Phase A: Preparation of ground base study plan 

Phase B: Laboratory test program 

Phase C: Experiment analysis program 

Phases A and B wete completed in 2*5 and 12 months, respectively, 
from the date of contract award. Phase C will be completed in 
19 months. 


SECTION II: DEVELOPMENT OF M.I.T. EFFORTS 


Turing Phases A, B and C the M.I.T. efforts basically followed 
the plan described in the proposal dated April 13, 1972 (Reference 1). 
The emphasis of the M.I.T. efforts was generally placed on the anal- 
ytical study. 

Phases A, B and C are described as follows: 

A. PHASE A: PREPARATION OF THE GROUND 

BASE STUDY PLAN 

The efforts during Phase A covered the following steps: 



J . - 


i 

Step 1: Collect currently available reports pertinent to this 

research 

} 

1 

Step 2: Consider all possible variables involved in the pro- 

posed research 

Step 3: Conduct rough, preliminary analysis , 

Step 4: Evaluate current experimental programs 

Step 5: Suggest modification of the current experimental pro- 

grams 

Step 6: Develop detailed plans of suggested experiment 

Step 7: Develop detailed plans of analytical research 

Efforts also were made to exchange information with researchers in 
other organizations assigned to studies other than thermal analysis. 

B. PHASE B: LABORATORY TEST PROGRAM 

Phase B was to conduct the study program defined in Phase A. 

Most of the M.I.T. efforts in Phase B were addressed to the heat 
flow analysis of the M551 experiment. 

A limited thermal analysis was also conducted by Teledyne Brown 
Engineering previously (Reference 2). The analysis used a finite 
difference computer program originally developed by Chrysler Corpor- 
ation. After that, other comprehensive thermal analyses were con- 
ducted by M.I.T. (References 3 and 4). Mathematical models based on 
the finite-element method were developed and applied to analyze the 
non-stationary heat flow of sample problems. Computers were used 
extensively in the analysis. 

The M.I.T. efforts in Phase B basically consist of the follow- 
ing steps: 

Step 8: Develop prelimxnary model 

Step 9: Improve analytical model 

Step 10: Analyze data obtained in ground base experiment 

In Step 8, a two-dimensional finite-element model was developed. 

The heat source was assumed to be concentrated at a point, to be dis- 
tributed evenly along the thickness direction, and to move along a 
circle at a constant speed. This analysis was expected to give 


2 



sufficiently accurate information about heat flow in most regions of 
the disc except the immediate vicinity of the heat source because 
the temperature distribution in regions apart from the heat source 
was presumably two-dimensional. However the temperature distribu- 
tion in regions near the heat source was predicted to be three- 
dimensional and the actual heat source was not concentrated at a 
point but distributed over some region of the disc surface. There- 
fore Step 9 was to develop a three-dimensional model and analyze the 
sample problem by use of the combination of the three-dimensional 
and two-dimensional models. 

The combination of those models yielded better information on 
heat flow in the vicinity of the heat source. 

Data obtained in a ground base experiment which was conducted 
by M.S.F.C. was sent to M.I.T., which was informed that the experi- 
ment was carried out under almost the same conditions as those des- 
cribed in Teledyne Brown Engineering's report (See Reference 2). 

The results obtained by the two-dimensional analysis did not show 
good agreement with the experimental results. As anticipated, the 
combination of two-dimensional and three-dimensional models showed 
the improvement of the analytical results, which were verified by 
the experimental results. 

C. PHASE C: EXPERIMENT ANALYSIS PROGRAM 

In experiments which have been conducted in space so far, a 
plan to obtain data on heat flow has not been included. Therefore 
the efforts during the remainder of Phase B and also Phase C have 
been addressed to improve the analysis and to study more details of 
the temperature distribution in regions near the heat source. 

Development of the analysis is presented in the following chap- 
ter. 


SECTION III; DEVELOPMENT OF ANALYSIS 


For the heat flow analysis of the M551 experiment, two types of 
mathematical models were basically developed: two-dimensional and 

three-dimensional models. 

At the beginning of Phase B, the two-dimensional model based on 
the finite element method was developed and applied to the heat flow 
analysis of the disc. Then, considering the varying thickness of 
the disc and the characteristics of experimental conditions, the 
three-dimensional model was developed. The model was combined with 
the two-dimensional model and used to obtain more reasonable results 
in the vicinity of a heat source. 


3 



In the following sections one of the typical sample specimens 
of the M551 experiment and its analytical results are described. 

A. SAMPLE SPECIMEN 

For the M551 experiment three types of materials were prepared. 
Their thicknesses were all different. The emphasis of M.I.T. analy- 
sis was placed on one of the materials — Stainless steel. The sample 
specimen used for the analysis is described as follows: 

1. It is a circular disc of 6.5 inches in diameter. 

2. It consists of four quadrants with different thicknesses 
of 0.020, 0.050, 0.125, and 0.250 inches. 

3. It has a back-up plate which shares the same shaft with it. 

4. It is stored in a chamber with other parts such as the back- 
up plate and a motor housing. 

The following values of the 321 Stainless steel are adopted as 
thermophysical properties of the specimen: 

Thermal conductivity 
Specific heat 
Density 
Melting point 
Latent heat of fusion 

The first three properties are the average values of the 321 stainless 
steel for the temperature range of interest. 

B. PREDICTED EXPERIMENT .L CONDITIONS 

Experimental conditions were assumed for the analysis as follows: 

1. A constant heat is generated by electron beam with 20 kilo- 
volts and 80 mllllamperes.* 

2. The efficiency of the electron beam is 80%. 

3. The electron beam diameter is 0.125 inch. 

4. The disc is heated along a 4.5-lnch diameter circle. 

5. The disc rotates clockwise with the relative velocity of 

37 Inches per minute against the electron beam through three quad- 
rants and then dwells for 30 seconds at the middle point on the weld 
circle in the 0.25- inch thick quadrant. 


= 12.9 Btu/hr/ft/°F 
= 0.14 Btu/lb/°F 
= 494 Ib/ft^ 

= 2600°F 
= 117 Btu/lb 


* From results of the ground-base experiment, 50 milllamperes 
were adopted in the analysis which used the combination of two- and 
three-dimensional elements. 
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TWO-DIMENSIONAL ANALYSIS 
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1. Two-dimensional Model . Figure 1 shows a mesh pattern and 
a node arrangement for the finite-element analysis of a disc. The 
model consists of two-dimensional finite elements. It contains 232 
elements and 216 nodes. 

In this analysis the following assumptions were adopted: 

a. Temperature is uniform in the thickness direction. 

b. Thermophysical properties are invaliant with temperature. 

c. The radiation between the disc and the surrounding parts Is 
neglected. 

d. The effect of solidification is neglected. ^ 

e. The Initial temperature of the disc is uniform at 0 F.* 

f. The heat source is concentrated at a point and moves along 
the weld circle. 

g. The effect of the metal melting on heat flow is taken into 
account by adjusting the intensity of the heat source. 

The following melting regions and intensities of heat source 
are assumed from results of the preliminary analysis. 


In Quadrant 1: 


1. The assumed melting region a band with a 0.475-inch 

width along the weld circle. 

2. The applied heating rate 1.017 Btu/sec. 

In Quadrant 2: 


1. The assumed melting 
width along the weld circle. 

2. The applied heating 


region a band with a 0.125-inch 

rate 1.084 Btu/sec. 


In Quadrant 3: 


1. Melting does not occur except at Node 113 in Figure 1. 

2. The applied heating rate 1.213 Btu/sec. 


* Basically, this assumption is permlssable because the linear 
analysis is adopted. In other words, the temperature, which starts 
from the initial temperature except zero, is obtained just by adding 
the initial temperature to the temperature calculated. 
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In Quadrant 4 : 

1. The assumed melting region a circle of 0.125-inch 

diameter (the size of the electron beam). 

2. The applied heating rate 1.110 Btu/sec (during a 

period when the region goes on melting) 1.213 Btu/sec (dur- 

ing the rest of the above period) . 

2. Analytical Results . Analytical results of the sample disc 
are presented in Figures 2 to 8. In those figures number*: on temper- 
ature history curves correspond with nodal names in Figure 1. The 
following results are shown in the figures: 

a. The temperature histories at points on the weld circle of 
4.5-inch diameter in Quadrants 1 to 3 (Figures 2 to 4) 

b. The temperature histories at the point where the electron 
beam dwells and at two other points in Quadrant 4 (Figure 5) 

c. The temperature histories of nodes located inside of the 
weld circle in Quadrants 1 to 3 (Figures 6 to 8) 

From those results the following items are pointed out: 

a. The melting occurs along the weld circle in Quadrants 1 
and 2. 

b. At the point where the electron be s dwells, the melting 
begins about 4 seconds after heat is supplit i there. 

c. The melting does no*, occur at points which are 0.325 inch- 
apart from the weld circle. 

d. The maximum temperature obtained by the two-dimfensional 
model exceeds 10,000 F in Quadrant 1. 

3. Discussion . The following discussions are probably made 
on the results described above: 

a. A point heat source assumed for the analysis is not realis- 
tic. 

b. Three-dimensional analysis should be applied to the region 
near the heat source in Quadrants 3 and 4. 

c. Heat radiation should be considered in the analysis. 

D. THREE-DIMENSIONAL ANALYSIS 

1. Three-dimensional Model . Three-dimensional finite elements 
were developed to idealize the disc as a three-dimensional model. 

In practice, the combination of two- and three-dimensional elements 
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were used to represent a mathematical model of the sample disc. The 
three-dimensional element was used for the region near the heat 
source and the two-dimensional element for the rest of the above 
region. 

A mesh pattern and typical nodal names used are shown In Fig- 
ure 9. The model contains 448 nodes and 336 elements. 

In the analysis the following assumptions are adopted: 

a. Temperature In the thickness direction is three-dimensional 
in the region of a band of a 0.5- inch width along the weld circle 
through Quadrants 2 to 4. 

b. In the rest of the above region temperature is uniform in 
the thickness direction. 

c. Thermophysical properties are invaliant with temperature.* 

d. The heat exchange due to the radiation between the disc and 
the surroundings Is not considered.* 

e. The effect of metal melting and solidification is not con- 
sidered. 

f. The initial temperature of the disc is uniform at 78 F. 

g. The heat source is uniformly distributed over the region of 
0.125-inch width band on the upper surface of the d .sc. 

2. Analytical Results . Results obtained by the three-dimen- 
sional analysis are shown in Figures 10 to 13. Temperature history 
curves in the figures denote the average of the temperatures at the 
nodes shown .in the parentheses. 

The following results are shown in the figures: 

a. The temperature histories at points on the weld circle in 
Quadrants i to 3 (Figures 10 to 12) 

b. The temperature histories at the point where the electron 
’am dwells and at two other points in Quadrant 4 (Figure 13) 

c. The temperature h.istorles of nodes located outside of the 
weld circle in Quadrants 1 to 3 (Figures 14 to 17) 

d. The temperature histories in the thickness direction in 
Quadrants 1 to 3 (Figures 12, 13, 16, and 17) 

The following items are observed from the figures: 

a. The temperature difference between the upper and the lower 
surfaces of the disc disappears in a few seconds after the heat 
source passes over. 


* The computer program is available for the an..iiysis consider- 
ing these effects. 
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b. The melting occurs along the weld circle in Quadrant 1 and 
at the point where the heat source dwells in Quadrant 4. 

c. The maximum temperature obtained is about 4,500 F in Quad- 
rant 1. 

3. Discussion . 

a. It is enough to use three-dimensional finite elements to a 
limited region close to the heat source in Quadrants 3 and 4. 

b. Even in the region in Quadrant 3 mentioned above, tempera- 
ture becomes uniform in the thickness direction of the disc in a few 
seconds after the heat source passes over, 

c. The temperature difference of the points close to the heat 
source in Quadrant 4 keeps almost constant value while ’leat is sup- 
plied. 

d. The heat radiation should be considered in the analysis as 
the next step because of considering the effect of the surroundings of 
of the disc. 


E. EXPERIMENTAL DATA AND ANALYSIS 

Data obtained by a gound-base experiment was sent to M.I.T. Ly 
Lockheed. It was informed that the experiment was basically con- 
ducted in the same conditions that were used in the Teledyne Brown 
Engineering's report (Reference 2). 

Generally speaking, the experimental results did not show good 
agreement with the analytical ones, as expected. This will probably 
come from two reasons, one of which is that details of the experi- 
mental conditions are not confirmed. The other reason is that fac- 
tors neglected or simplified in the analysis have influence, which 
is not predictable, on the analytical results. 

Under the circumstances three nodes cJose to the weld circle are 
chosen to compare the experimental results with the analytical re- 
sults. 

Figure 18 shows the temperature history of a point in Quadrant 1. 

The experimental result Is different from the analytical one. 

The difference is possibly caused by the differences between the heat 
Intensities supplied in the experiment and those assumed for the 
analysis. In addition, the effect of melting on the heat flow of the 
disc will not be negligible because the thickness of the disc Is very 
thin in this quadrant. The velocity of the electron beam is slower 
than that adopted in the analysis. 

Figures 19 and 20 show the temperature histories of typical 
points in Quadrants 2 and 3, respectively. In those figures, the 
correlation between the experimental and the analytical results 
seems to be reasonable. However the heat actually supplied to the 
disc is presumably more than that used in the analysis. The effect 
of the heat radiation is also observed from Figures 19 and 20. 
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F. COMPUTER PROGRAM 

Computer programs based on the finite element method were de- 
veloped for the analysis of non-stationary heat transfer problems. 
Among various programs two programs were mainly used to analyze heat 
flow of the specimen. One Is for the two-dimensional analysis and 
another for the three-dimensional analysis. In the three-dimensional 
analysis program a two-dimensional finite element as well as a three- 
dimensional one are prepared to Idealize the specimen as the combina- 
tion of the two types of element. 

The following guidelines were considered in developing the 
computer programs. 

1. The program can be applied to a heat flow analysis of any 
two- or three-dimensional body. 

2. Both specified temperature and supplied heat are treated as 
boundry conditions. 

3. The effect of radiation can be considered. 

4. Metal melting and solidification can be considered. 

5. Non-stationary heat flow problem as well as stationary one 
can be solved by the s„me program. 

6. Thermophysical properties are allowed to depend on tempera- 
ture. 

7. The program should be flexible for future improvement or 
modification. 

These Items have been achieved so far. Details and manuals of the 
programs will be presented later. 

G. CONCLUSION 

The heat flow analysis of the M351 experiment wes conducted at 
M.I.T. In the analysis, only the disc was considered. In other 
words, the surroundings of the disc were ignored. The computer 
programs based on the finite-element method were developed for the 
analysis. 

The disc was idealized as the two-dimensional mathematical 
model at first. Then the three-dimensional model was developed, 
which consisted of the combination of two- and three-dimensional 
finite elements. 

The analytical results obtained by the three-dimensional model 
were compared with the limited experimental data. From the compar- 
ison it became clear that the heat radiation had influence on the 
heat flow of the problem in consideration. 

The computer programs are now available for the analysis con- 
sidering the radiation. The programs also cover the effect of metal 
melting and solidification and the temperature dependency of thermo- 
physical properties. 
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SECTION IV: APPENDICES 


APPENDIX A 


Symbols 



X. yt 2 


£ 

\ 

P 

O 


Area 

Heat Capacity 

Absolute temperature 

Thermal conductivity 

Components of unit normal vector 

Rate of internal heat generation 

Temperature 

Prescribed temperature 

Melting point 

Time 

Volume 

Position coordinates 
Emlsslvlty 

Latent heat of fusion 
Density 

Stefan-Bolzmann Constant 



APPENDIX B 


Generally speaking, the finite element method Is based on 
the principle of virtual work or the principle of complementary 
virtual work. They are dual principles. 

In this report, the principle of virtual work is applied to 
the analysis of a non-statlonary heat transfer problem. This 
principle is widely applied to analyze problems in various fields 
because the existence of potential function is not required. 

For the analysis of the MSSl experiment, both two- and three- 
dimensional models were developed by use of Che finite element 
method. Emphasis of the formulation is placed on the three-dimen- 
sional model In Che following section. 

The governing equation and the corresponding boundary conditions 
are given for the three-dimensional heat transfer problem as 
follows: 


_ 3 . 

3x 


k ^ 1 

+ ± 


+ ± 

Ik 51 1 

[x3x J 

3y 

ly3y 

3z 

rz3z J 




in V 


(A-1) 


nk l^+nk |~+nk ~ + ea(K^ - K^) 
X X 3x y y 3y z z 3z ' o 


on 


on 


where A^ denotes the surface on which temperature is 


A 2 (A-2) 

A 3 (A-3) 

specified 


and A 3 the surface from which heat emission due to radiation occurs. 

Eqxiatlon (A-3) presents the Stefan-Boltzmann law modified for a 
grey body. The principle of virtual work which is equivalent to 
equations (A-1) and (A-3) is obtained as follows: 



3T 


3z 


+ Cp|^ 6 t| dV 


(A-4) 


In equation (A-A), temperature must satisfy equation (A-2) on 
boundary A 2 « In other words, equation (A-2) is the subsidiary 

condition for equation (A-4). 

In addition, it is noticed here that terms concerning time such 
as ^ and * are assumed to be constant during variation 6 t in deriva- 
3t " 

tlon of equation (A-4) . 
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APPENDIX C 


Finite Element Approach 

In formulating a non-statlonary heat transfer problem, the 
concept of finite element Is adopted for the space coordinate 
and that of finite difference for the time coordinate. 

As seen In equation (A-4), this equation leads to the nonlinear 
one. To solve the equation, an Iterative procedure Is commonly 
used. However, If that procedure Is applied to a transient problem 
such as the non- stationary heat transfer problem, computing time for 
the analysis will Increase tremendously. In order to avoid this 
disadvantage and also to keep the simplicity of the finite element 
method In Its formulation, temperature which determines the emission 
of radiation energy from a body Is assumed to be known. In other 
words, the emission Is assumed to be obtained by using the tempera- 
ture at the latest time step which Is known. This assumption will 
be valid If the temperature change during one time step Is small. 
Under the above assumption, the formulation of the problem by 
use of the finite element method becoises straightforward. 

As shown In Figure A-1, the two-dimensional and three-dimension- 
al elements are used together to analyze the heat conduction problem 
of the disc. However, three-dimensional analysis Is mainly described 
here. 

In the finite element analysis, a body of Interest Is Imaglnarl- 
ly divided Into a finite number of elements. Therefore, equation 
(A-4) Is written as follows: 




+ k ^ 6 

y 3y 



+ k 


z 3z 


6 


3z 

V / 


I 




”2 

+ 2 //. £a(K^ - K^) (STdA 

1=1 ^ 2 ^ ° 


N 


- ^ ///v ^ 


(A-5) 


1-1 


where N and N^ are total number of elements and that of areas con- 
cerning radiation; respectively. 

For a typical three-dimensional element shown In Figure A-2, 
temperature dlstribotlon In It is to be assumed by the followiug 
polynomials . 
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(A-6) 


T ■ a + tx.x + a*y + gl.z 
0 12 J 

This is also written in the matrix form as follows. 


[1, X, y, z] 


a. 


a. 


a. 


Equation (A~7) is further written concisely 
T - [A] {a} 


(A-7) 


(A-8) 


where 

[A] - [1, X, y, z] (A-9,1) 

{a}^» [Uq, a^, a^, a^] (A-9.2) 


As seen from Figure A-2, {a} is related with the temperatures at 
the nodes of an element as follows. 
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^1 
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“l 

^3 

1 s 

1 

*3 

^3 

^3 

1 

“2 



1 


^4 



“3 









. / 


(A-10) 


where Xj^, y^^, .... and z^ denote the coordinates of the corres- 
ponding nodes. 

For simplicity, equation (A-10) is written as follows. 
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{T} * [N]{a} 


(A-11) 


where 

* [T^, T^. T^, T^] 


(A-12.1) 


[N] = 


1 Vi 


^ ^2 ^2 "2 


^ ^3 ^3 "3 


1 ^ V4 "4 


(A-12.2) 


Substituting equation (A-11) into equation (A-8) , the assumed tem- 
perature distribution in an element is represented in the following 
form. 


T * [A] [N]“^{T) 


(A-13) 


By using equation (A-13) , each Integral term in equation (A-5) is 
easily written in terms of nodal temperature. 

Differentiating equation (A-13), the following expression is 
obtained. 


9x 


3T 

3y 


[B] [N]“^{T) 


3z 


(A-14) 


where 
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0 10 0 


[B] 


0 0 10 


0 0 0 1 


(A-15) 


Using equation (A-14) , the first three terms in equation (A-5) are 
written as follows. 


6{T}^ ([N]"V///y [a]^[K] [B]dV[N]~^{T} 

i 

where 



(A-16) 


(A-17) 


Similarly, the fourth term is written as follows. 

/// Cp[A]^[A]dV [N]“^{T> (A-18) 

where 





3T, 3T, 

3, 4 


>c 

9t at 


(A-19) 


The fifth term in equation (A-5) reminds us of the assumption men- 
tioned above. Under the assumption that the absolute temperature 
on the boundary A„ is knovm, that the term is represented as 

^ j 


follows. 
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f 


6{T}^ UNf^)'’^ //^ tAl^ea(K“ - K^)dA 


(A-20) 


The last term Is written In a similar form to that of equation (A-20) . 


6{Tr ([N]"^)"^ JjL [A]"QdV 

1 


.t: 


(A-21) 


For the expressions (A-16) , (A-18) , (A-20) and (A-21) , the following 
matrices and vectors are defined. 


[H] - ([N]~S^ JJL [B]"tK][B]dV[N] 


-1 


(A-22) 


IP] - ([N]"^)*^ IIJv^ Cp[A]^ [A]dV[N] 


-1 


(A-23) 


« ([N]"^)”^ JJa^ lA]^ eo(K‘* - K^)dA 


iT /••A , rA \ 


(A-24) 


(F } = ( m "^)^ ///v^ lA]^ QdV 


T 1 


(A-25) 


Using the expressions (A-22) to (A-25) , equation (A-5) is written as 
follows. 

N . N T ”2 

Z 6{T> [H]{t} + I alTl’^lPHT} + Z 6{t} {F } 
l-l 1-1 1-1 


- Z ip } - 0 (A-26) 

.■1 " 
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The summation symbol In the above equation denotes that components 
In a matrix or vector concerning each element are entered In their 
position of the global matrix or vector. For simplicity, the same 
symbols as those for an element are adopted to represent the ch;>r- 
acterlstlc matrices and-vectors for a whole body. Then, equation 
(A-26) Is written as follows. 


6{t}^[H]{t} + 6{t}'^[P]{t} + 6{t}'^{f_} - 6{t}^{F-} - 0 

(A-27) 

In equation (A-27),6{T} Is arbitrary. Therefore, the following 
equation Is obtained so that equation (A-27) Is held, 


[H]{T} + [P]{T} - {F} - {0} 


(A-28) 


where 


{f} - {Fq} - iFj^} (A-29) 

To the time coordinate, a finite difference method Is usually applied. 
In this report, the following finite difference scheme Is used. 


{T>i - ^ ({T}^ + (T}^_^) 


(A-30) 


where subscripts 1 and 1-1 denote time, t^, and t^ respectively; 
At Is the time interval between t^ and t^ Using equation (A-30), 
equation (A-28) Is transformed as follows. 

([H]+ ^ [P])A{t}^ - -2 [H]{t}^_^ + + {f}^ (A-31) 


where 



I 


i 

! 

t 


A{T}^ = {t}^ - 


(A-32) 


Therefore, if the Initial condition is given, equation (A-31) can 
be solved step-by-step at every time increment. 






APPENDIX D 


Finite Element Models 

Concerning the nest conduction analysis of the disc, several 
types of finite element were used, as shown In Figure A-1. Two of 
them are two-dimensional models and others are three-dimensional 
models. Both of the two- and three-dimensional models contain a 
composite element. These composite elements are introduced in order 
to make it easy to generate data input for the analysis and also to 
Improve the accuracy of the solution. 

As seen from Figure A-1, the two-dimensional composite element 
consists of the average of four triangular elements and the three- 
dimensional element, the average of ten tetrahedra. In addition, 
these composite elements as well as a simple triangular element or 
a simple tetrahedron assure the continuity of the assumed temperature 
distribution on the boundary between adjacent elements. This contin- 
uity condition is required so that the solution of equation (A-5) 
approaches the exact solution when the number of elements increases. 



APPENDIX E 


A Basic Concept to Consider Metal Melting and Solidification 

In the finite-element analysis of heat conduction problems, It 
Is possible to consider metal melting and solidification. For exam- 
ple, metal melting is considered In the analysis »s follows. The 
problem Is to estimate how much heat is spent for metal melting and 
what Is the temperature of an element of interest. 

First, consider an element with the volume of dV. Tb^r. , suppose 
the average temperature of the element is T, at time id T is 

below the melting point of the metal. If the tempers . Increases 
during the next time Increment, we may expect one of ‘ * cases as 

shown in Figure A-3. 

In Case a) , the average temperature of the element is still 
below the melting point. Therefore, the same procedure as before 
will be taken for the following calculation. 

In Case b) , the average temperature Is beyond the melting point 
and heat storage during the time Increment Is enough to melt the 
element completely. Therefore, we ha-'e the following relation. 


‘‘t “ “l ^2 

where and denote heat stored in the element and heat spent to 
melt the clement, respectively, q^^ and q^ are represented as follows. 


q^ - Cp (T^ - Tj^)dVAt (A-34.1) 

q^ - Cp (T^ - T^)dVAt (A-34.2) 


In equation (A-33), q* is calculated from the temperature, T^, which 
Is the original one obtained at time, t 2 > as follows 

q^ - Cp(T^ - T^)dVAt (A-35) 

Using the latent heat of fusion. A,, q„ is obtained as follows 

t M 

q„ - dVA^ (A-36) 
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Therefore, substituting equatlous (A-34) to (A-36) into equation 
(A-33) , the modified temperature, at the time t 2 can be calcu- 
lated. For the next calculation, the temperature T 2 is treated as 
the initial condition. 

In Case c) , the average temperature of the element interest 
is also beyond the melting point, but heat storage during the given 
time increment is not enough to melt the whole element. Therefore, 
the following relation is obtained instead of equation (A-33) 


and q^ are given in the same form as equa-ions (A-35) and 
(A-34.1), respectively. On the other hand, q^^ is given as follows. 

q^ - dV'A^ (A-38) 

where dV' < dV which means that the whole element of interest does 
not melt. For the next calculation, the melting point T^^ is treate 

as the initial condition. 

In the case of solidification, the procedure considered in the 
cjmputer program Is the same as that for melting, except that the 
temperature decreases. 
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Figure 8 Calculated Temperature Histories of Nodes inside of Weld Circle in Quadrant 
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FIGURE 9 Mesh Pattern and Node Arrangement 

for F...tlte Element Analysis of Disc 
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Figure 10 Calculated Temperature Histories of Points on Weld Circle in Quadrant 





u 

UJ 

w 

z 

o 

I- 

< 

H 

z 

H 

Z 

UJ 

S 

cr 

UJ 

Q. 

X 

UJ 

2 

o 

ir 

u. 

UJ 


I- 


(Jo) 3ynivy3dW3i 





0) 



VI 



d 



4J 



cO 


CN 

M 



0) 


44 



a 

a 



0) 


u 

44 


T> 



Cd 



3- 

4i 





B 



■H 

o 





Qt 

ta 


iH 



U 

"9 


Vi 

n 


•M 

3 


O 

« 


•3 

at 


iH 

at 


at 

at 


;* 

JB 



tj 


B 

0 


o 

at 



h 


te 

a 


v> 

S. 


B 



•H 

la 


O 

Tt 


a 

« 


4-1 

•3 


o 

o 


po 

B 


at 



■H 

<U 


Vi 

o 


0 



Vi 

a 


K 

at 


•rt 

a 


X 

«* 


at 

at 


Vi 

3 

s. 


U 

4 

& 


Vi 

w 


at 



i* 

at 


§ 

60 


H 

at 



lx 


•3 

01 

> 

a 

> 

ot 

u 

Id 

u 

cd 


cd 

I — 1 


MX 

3 

Ja 

Vi 

U 

♦J 

3 

rH 


0) 

at 



o 

u 

Vi 


Q 

(U 


5 

a 


01 

a 

rH 

T) 

3 

rH 

ID 

01 


01 

M 

3 

> 

4J 

p 

U 


o 

a 

44 

M 

u 

O 

a 




at 



«j 



o 



a 




32 


I 

Ii4 

IS 



(do) 3dnivy3diN3i 


0 
u 

0) 

u 

1 

s 


<*> 



(B 

« 6 

^ t 

Si 


U O 


SI ■ 

^ S) 
0 M 
h 0 
0 4i 
> 0 
0 

0 

0 ‘ 
s: 


0 » 
tj *j 
0 

0 .A 
0 

t § 

S 0 

L> T* 


0 

U 

o 

z 


FIGURE 12 Calculated Teatperature Histories of Points on Weld Circle In Quadrant 
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Figure 14 Calculated Temperature Histories of Nodes outside of Weld Circle in Quadtant 
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FIGURE 16 Calculated Tesperature Histories of Nodes outside of Weld Circle in Quadrant 
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FIGURE 17 Calculated Temperature Hletoriea of Nodes concerning Dwelling Heat 



















(a) BASIC TRIANGULAR ELEMENT 




(c) BASIC TETRAHEDRON 



(d) COMPOSITE ELEMENT 


Figure A-1 Finite Elements Prepared for Analysis 
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Note: (a) Heat is not supplied enough to melt an element. 

(b) Heat is supplied enough to melt an element. 

(c) Heat is supplied enough to melt part of an element. 


Figure A-3 Models to Consider Metal Melting 
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